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The presence of sound waves in one or the other of the fluid regions 
on either side of a shock wave is made apparent, in the region under 
superpressure, by acoustic waves (reflected or refracted according to 
whether the incident waves lie in the region of superpressure or of sub- 
pressure) and by thermal waves. The characteristics of these waves are 
calculated for a plane, .progressive, and uniform incident wave. In the 
case of refraction, the refracted acoustic wave can, according to the 
incidence, be plane, progressive, and uniform or take the form of an 
"accompanying wave" which remains attached to the front of the shock 
while sliding parallel to it. In all cases, geometrical constructions 
permit determination of the kinematic characteristics of the reflected 
or refracted acoustic waves . The dynamic relationships show that the 
j amplitude of the reflected wave is always less than that of the incident 
wave. The amplitude of the refracted wave, whatever its type, may in 
certain cases be greater than that of the incident wave. 


1. BASIC CONSIDERATIONS 


1.0. Generalities 

1.0.0 Subject of the report .- The reflection and the refraction of 
sound waves by a shock wave are of interest not only in acoustics. These 
phenomena may also be utilized in aerodynamics to detect or measure the 
development of shocks. 

The present paper tests the theory by application to a relatively 
simple case: that of a plane and uniform shock wave separating two 

regions in both of which the fluid is in pressure and temperature equi- 
librium, except for the acoustic phenomena) the incident acoustic wave 
is plane, progressive, and uniform. 

"Reflexion et refraction d'ondes acoustiques par une onde de choc. 
Acustica, vol. 5, no. 1955? pp. 1^9-165* 


It 
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This problem has already been treated in several publications. 

Burgers (ref. l) and Blokhintzev (ref. 2), on the one hand, have studied 
analytically the case of normal • incidence in a fluid without viscosity. 
They show that, in order to satisfy the boundary conditions, it is 
necessary to introduce, in addition to the reflected or refracted acoustic 
wave, a so-called "entropy" wave. This nonprogressive wave involves tem- 
perature fluctuations. It is the asymptotic form of a thermal wave in 
the expression of which the damping terms have been eliminated. The 
results regarding this wave thus obtained apply therefore only in the 
neighborhood of the shock wave. 

Sauer (ref. 3), on the other hand, indicated the conclusions to be 
drawn from Huyghens ' construction and demonstrated the existence of 
incidences for which no refracted wave of a usual type (that is, a purely 
longitudinal wave) can exist. 

The purpose of the present report is to examine the points left 
undecided by the preceding authors, namely: the behavior of the thermal 

wave at large distance, and the nature of the phenomenon if there is no 
refracted wave of longitudinal type present. We have thus been led to 
treat the problem in a more general form: we made an analytic study, 
taking into account the viscosity and the caloric conductivity, for 
arbitrary incidence. 

1.0.1 Basic hypotheses and approximations .- (a) In the calculation^ 
we shall consider the shock wave as a discontinuity. The thickness and 
the duration of the development of the shock are actually very small 
compared to the wave lengths and the lengths of the acoustic periods. 

(b) The equations governing the shock wave in a steady state remain 
valid in the presence of acoustic waves. The slowness and the small 
amplitude of the acoustic phenomena compared to the shock phenomena 
justify this hypothesis, at least in the first approximation. 

(c) In the expression of the boundary conditions, we neglect the 
small oscillations of the front of the shock wave about its mean position. 
This is an approximation which is justified- in all acoustic problems. 

Due to the fact that the .acoustic displacements are very small with 
respect to the wave lengths, the errors thus committed are of the second 
order . 


^This hypothesis is implicitly contained in Burgers calculation. 
This author calculates the motion of the shock wave by a method which is, 
incidentally, not applicable to the case of oblique waves; then, elimi- 
nating in the boundary conditions the terms of the second order, he 
achieves a simplification which corresponds to our hypothesis. 
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The "boundary conditions will then he expressed in point form, that 
is, that the relationships expressing them contain only the values of 
the quantities concerned taken from hoth sides of the boundary (front 
of the shock wave) at one point of that boundary. 

1.0. 2 Irreversibility of the shock wave and its consequences .- The 
development of a shock wave is irreversible. Thus it is evident before- 
hand that it will not be possible to satisfy the boundary conditions if 
in the two regions separated by the shock wave only acoustic phenomena, 
that is, adiabatic phenomena, are introduced. 

As a result, the small motions produced by the meeting between a 
shock wave and an acoustic wave will be expressed in their general form 
as it results from the work of Stokes and Kirchhoff, that is, taking 
into account the viscosity and the thermal conductivity of the gases. 

1.0. 3 Separation of the boundary conditions .- Since these conditions 
are in point form, one can divide them into two groups. 

(a) Kinematic conditions which govern the nature and form of the 

waves produced by reflection or refraction: When all the waves concerned 

are plane, progressive, and uniform, these conditions may be put in tfce 
geometric form of Huyghens 1 construction. In the other cases it is 
p.eces 8 ary to resort to calculation in order to determine the nature and 

^the form of the waves produced. This is what is done for instance in 
optics in order to treat the problem of total reflection and to obtain 
the Cauchy wave. 

The kinematic conditions simply express the fact that on the 
boundary surface, the development in space and time of the phenomena 
associated with the reflection or the refraction is the same. Thus it 
is not necessary to use, for expressing them, the theoretical relation- 
ships which govern the shock wave. 

(b) Dynamic conditions which come into consideration only for 
calculating the amplitude of the phenomena whose nature has been deter- 
mined by the kinematic conditions. 

We shall study these two groups of conditions separately. 


1.1. Characteristics of the Shock Wave. 

Boundary Conditions 

1.1.0 Definitions, notation .- In designating the two regions of the 
fluid separated by the shock wave we shall avoid using the terms upstream 
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and downstream which. In aerodynamics as well as in hydraulics, give rise 
to ambiguity. If the shock wave -is motionless (wind tunnel, hydraulic 
jump),- the region downstream is under superpressure. If the wave is 
propagated in a fluid at rest (shock tube, tidal wave), the upstream 
region is under superpressure. 

Thus our notation will be 

Eq = region under subpressure, extending from S toward the positive x 

= region under superpressure, extending from S toward the 
negative x 

S = front (plane) of the shock wave separating the two regions, per- 
pendicular to the axis of the x 

We shall use three systems of axes: 

0XQyz fixed with respect to the fluid in the region Eq 

ox^yz fixed with respect to the fluid in the region Ej_ 

oXyz fixed with respect to the front of the shock wave S 

With 

&q > 0 being the velocity of propagation of the front S in the fluid 
contained in Eq 

> 0 the velocity of propagation of S in the fluid contained in E-j_ 

U the difference between the flow velocity of the fluid in E-j_ 

and its flow velocity in Eo 

We have: 

Xi = Xq - Ut X = Xq - aQt = x-j_ - ajt (l.l.l) 


and 


&0 ” = u 


( 1 . 1 . 2 ) 
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Finally we define: 


Region Eq Region E-j_ 


Pressure 
Specific mss 

Absolute temperature 
Velocity of sound 


P0 

*0 

c 0 


(1.1.3) 




We shall take as the independent variable defining the amplitude 
of the shock wave the ratio: 


i - *oh 


(1.1.4) 


' | decreases from 1 to 0 in proportion as the amplitude of the shock 
increases . 

We shall use dimensionless functions for describing the phenomenon. 
We give them below with their expressions as functions of | derived 
from the classical theory of shock waves. 

Putting 


H = 


_ 7 + 1 


7 - 1 


c 0 _ 

(*0 

f /2 


iln + il 

C 1 

V T i> 

/ 


1 + 

R - 

p o 

_ i 

+ n£ _ 1 

p i 


+ | p2 

; _U_ 


(m- ■ 

• DU - i) 


1/2 


(1.1.5) 


( 1 . 1 . 6 ) 


(1.1.7) 
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Aq = — - 
"0 c 0 


h + t 


11/2 


UTTTTI 


Aq> i 


( 1 . 1 . 8 ) 




1/2 


A]_ < 1 


(1.1.9) 


Finally, we shall have to use the following functions in which 
appear the derivatives R' and Y' of R and Y with respect to g: 


|Y' .= 


(h - 1 

)t + (H + 2)|_ 


2(|i + |) J 

(h + l)|(|i + |) 

]l/2 


(1.1.10) 


( y. - lWn + 1 R_|_ , 

IH + 1/ W - 1 S R 


fa - 1\ 2 / (n + i) 2 i 

i H + 1 j 1(1 + j-t i. ) ( M- + I) 


- 1 


(1.1.11) 


1_ 

Rr 


h + £ 


l(i + nl) 


[ 


1/2 


( 1 . 1 . 12 ) 


G = - 1 JIY* 

H - 1 


H + (h + 2)|Ax + 1 W 2 

2((i +|) \1 + n|/ 


(1.1.15) 
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H = 





n + (14 + 2)6 

O I “ t 

a + 1 

H + 1 

iJL + 1 

luTTT) 


1/2 


(1.1.14) 


Tables I and II give the values of these various functions, most 
frequently calculated according to rule which gives sufficient accuracy 
for our present purpose. 

1.1.1 Wave magnitudes .- The presence of waves with small amplitude 
which have an acoustic origin causes in the regions Eq and E^ small 
variations of the pressure, the specific mass, and the temperature about 
their mean values. 

In order to represent them we set: 

0 ) = 5P/P = relative excess of pressure 

S = Sp/p = relative excess of specific mass, or condensation (1.1.15) 

9 = 6T/T = relative excess of temperature 

In the case of usual' acoustic waves, d5, s, and 9 are very small 

and may be considered in the calculation as infinitesimals of the first 
order . 

With regard to the vibratory velocities, we shall have to consider 
only the component u along ox. 

1.1.2 Boundary conditions .- According to our hypotheses, the boundary 
conditions must express that the relationships (1.1.4) to (1.1.7) which 
govern the shock wave are maintained in the presence of oscillations of 
acoustic origin. 

By differentiation of (1.1.4), (1.1.5)* and (1.1.6) we obtain 


cE-l - cdq = -H/i 


(1.1.16) 
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.01 - B q = -2 £1 B4 = 2| - <5q) (1.1.17) 


S 1 - 8 0 - " r- 8 * - i - <5 0 ) (1.1.18) 

In order to obtain* starting from (1.1.7)* a correct expression 
one must take some precautions. One could be tempted - as we were - to 
assume* as it is usual in acoustics that the presence of the acoustic 
waves does not modify the sonic velocities Cq and c^ 

This approximation remains correct in the interior of the regions 
Eq and Ei but at their boundaries it c ann ot be assumedj the relation- 
ship (1.I.5) contradicts it. Thus one must differentiate (1.1.7) with 
the assumption that Cq is variable. One then obtains 

SU = e 0 8X + Y8c 0 (1.1.19) 

with 6U = Ui - Uq* the difference between the components Ui and Uq 
of the vibratory velocities being given by 

5Y = Y'5g = -gY'^i - a> Q ) 

On the other hand* cq is given by the classical relationship 

c o = 7 P o/ p o 

whence 

25c 0 /c 0 = 5P 0 /P 0 - BPo/Pq - - s 0 

Ab we shall show later on in our problem* only incident acoustic 
waves can exist in the region Eq. The oscillation phenomenon is 

therefore adiabatic which permits us to write s Q = a wV* hence 


5c o/ c o = + ^ 
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The relationship ( 1.1. 19) may then he written 


]i! tr 2 “ (ir+T + * Y ’H - 


( 1 . 1 . 20 ) 


It will he convenient to use this equation in the form 


=1 


0 - 


c l u l " G — “1 = pl c Q u 0 " H — 


( 1 . 1 . 21 ) 


where r, G, H are the functions defined, respectively, hy the 
relationships (1.1.5), (1.1.13), (1.1. 1*0. 


1.2. Waves of Small Amplitudes in Gases 

1.2.0 Mathematical representation . - The mathematical representation 
of this problem is a result of the work of Stokes and Kirchhoff. An 
outline for it can he found in Lord Rayleigh's hook (Theory of Sound, 
vol. II., section 247). We recapitulate it with our notation. 


u, v, 

P 

s = 5p/p 
P 

m = SP/P 
T 

0 = 5T/T 

7 

V 

k 


w 


components of the velocity of a particle 
specific mass of the gas at rest 
condensation 

pressure of the gas at rest 

relative excess of pressure 

absolute temperature of the gas at rest 

relative excess of temperature 

ratio of the specific heats 

coefficient of kinematic viscosity 

coefficient of thermal conductivity, accord ing to Kirchhoff 
equal to 5v/2 


sonic velocity 



10 


NACA TM 1409 


Equation of continuity 


ds/dt + du/dx + dv/dy + dw/dz =0 ( 1 . 2 . 1 ) 


Dynamic equations 

~ + ~~ - v(a^u - — ■) and similarly in y,z ( 1 . 2 . 2 ) 

at 7 ox \ 3 ot ox/ 

The logarithmic differentiation of the equation P = prT yields 
the equation of state 


(D = S + 0 (1.2.3) 

Since the motions are very small, the terms of the second order are 
neglected. The total relative temperature excess 9 is the sum of the 
excess 9 a , resulting from the adiabatic compression, and the excess 9^, 
stemming from the conductivity. Thus we have 


9 a = (7 - l)s d9-fc/dt = kA 2 9 

Hence we obtain finally the thermal equation 

d9/dt = (7 - l)ds/dt + hA 2 9 (1.2. 4) 

The calculation then is performed as follows: We eliminate co 

utilizing (1.2.3)- We then assume that the different variables depend 
on the time only through the factor e^ (h can be complex). The 
factor h therefore replaces the sign of differentiation with respect 
to time. 

Equations (1.2.1), (1.2.2), and (1.2.4) then are written 


du/dx + dv/dy + dw/dz + hs = 0 

( 1 . 2 . 1 ') 

vA^u = -c)M/cSx, and similarly in y,z 

( 1 . 2 . 2 ' ) 

h0 - kA 2 9 = (7 - l)hs 

(1.2.4’) 
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where we have set 


M = (c 2 /V + hv/3) s + (c 2 /r)^ (1.2.5) 

We eliminate s "between (1.2.1*) and (1.2.4*) which giveB 

du/dx + 3v/dy + dw/dz = -(h9 - kA 2 9)/(7 - D (1.2.6) 

nnrl M is written 

(1.2.5’) 


Let us differentiate equations (1.2.2* ), respectively, with respect 
to x, y, z j and add term by term. We take (1.2.6) into account and 
eliminate M, using (1.2.5' )j there remains 


M = 


- S~ 


r - 1 




hr 1 



(1.2.7) 


a particular solution of which is given "by the form 

9 = + ® 2^2 ( 1 . 2 . 8 ) 

where B-^, Bq are constants, and and Q 2 functions of x, y, z 

which satisfy the relationships 

A 2 Q]_ = A^2 = m 2^2 (l*2.9) 

in which xn^ and m| are determined by the equation 



( 1 . 2 . 10 ) 



12 


MCA TM 1409 


with a being the dimensionless number 


0 = hk/c 2 (1.2.11) 

The components u, v, w of the vibratory velocity then derive, 
except for the common factor e* 1 ^, from the potential Qi + 0.2’ ® Le 

constants and are taken from (1.2.6) which gives 

2 

B i = (7 - 1) ~ - (i = 1,2) (1.2.12) 

kmf - h — 

i 

The condensation s^ corresponding to Q,^ is 

Si = -(m 2 ^h)cpie ht (1.2.13) 

The relative excess of pressure d3j_ is obtained by ^combination of 
(1.2.3) and (1.2.13) 


co-r 


m 2 km| - ht 


(1.2.14) 


1.2.1 Simplifications, separation of the wave types (acoustic and 
thermal ) . - In the case where the origin of the waves is acoustic, we 
shall have h = lo - 6. In the relationship (1.2.10) which determines 
the ratios m/h, h appears only in the dimensionless number a = kh/c 2 . 

When this number is very small, the roots of (1.2.10) assume very simple 
approximate values. In air, in cgs units, v = 0.l46j thus k = 5^/2 = 0.3&5 
At the usual temperatures c is of the order of 34 X 10^ cm/s. Finally, 
if we restrict ourselves to the upper limit of an - already very ultran- 
sonic - frequency of 80,000 cycles per second, the circular frequencies to 

are smaller than 0.5 X 10^s - -*-. Under these conditions, kto/c 2 is less 
than 1.6 X 10“^. 

We are therefore justified in assuming that the number a will 
remain very small as long as the damping factor S is not very large. 

We shall therefore make the corresponding simplifications, except 
for verifying - when the occasion arises - that the condition a « 1 is 
always satisfied. 
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The roots of (1.2.10) assume then the approximate values 

m^ » (h/c) 2 « 7h/k (1.2.15), ( 1 . 2 . 16 ) 

The first gives us a velocity potential for an acoustic wave in 
its usual form. The second gives us the potential for waves of thermal 
conductivity . 

In the case of plane waves, with the coordinate z eliminated hy 
a suitable choice of axes, we shall be able to write 


inuCax+Py) 

% - e 1 

with a , 2 + 0 2 = 1 , in order to satisfy ( 1 . 2 . 9 )* 

We then obtain the table of the following values: 



Acoustic wave 

Thermal vave 

Q 

e +h(ox+ 0 y)/c 

e i\J 7 h/k( ax+py ) 

s 

-(h/c 2 )Qe ht 

-( 7 /k)Qe ht 

e 

-(7 - l)s 

-s 

03 

7 s 

0 

11 

teas 

+ \Jkh/ 7 as 

V 

±c(3s 

+\jkh/ 7 ps 


(1.2.17) 


( 1 . 2 . 18 ) 
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1.5* Numerical Tables of the Different Functions 
Introduced in the Calculation 


Table I 


5 

r 

R 

R' 

Y 

-|Y' 

*0 

A 1 

1 

1 

1 

0.714 

0 

0.714 

1 

1 

Q.9 

0.985 

0.95 

0.735 

0.0756 

0.726 

1.046 

0.956 

0.8 

0.97 

0.85 

O.76 

0.160 

0.737 

1.10 

0.910 

0.7 

0.95 

0.78 

0.78 

0.26 

0.756 

1.17 

0.86 

0.6 

0.93 

0.70 

0.80 

0.39 

0.78 

1.25 

0.81 

0-5 

0.90 

0.62 

0.83 

0.53 

0.81 

1.36 

0.76 

0.4 

0.87 

0.53 

0.85 

0.72 

0.85 

1.51 

0.70 

0-5 

0.82 

0.45 

0.88 

O.96 

0.92 

1.73 

0.63 

0.2 

0.75 

0.38 

O.91 

1.38 

1.04 

2.10 

0.56 

0.1 

0.62 

0.26 

0.94 

2.18 

1.35 

2.95 

0.48 

0.05 

0.48 

0.214 

0.956 

>27 

1.82 

4.15 

0.431 

0.025 

0.36 

0.191 

0.964 

4.75 

3.18 

5.86 

0.405 

0 

0 

0.1667 

l 

0.973 

OO 

CO 

00 

0.378 


Table II 


1 

J 

i/Rr 

H 

G 

1 

0 

1 

1 

1 

0-9 

0.0007 

1.094 

1.00035 

1.00035 

0.8 

0.0031 

1.21 

1.00168 

I.OOI65 

0-7 

0.0079 

1.36 

1.0052 

1.0044 

0.6 

0.016 

1-55 

0.012 

1.009 

0.5 

0.029 

1.80 

1.025 

1.018 

0.4 

0.051 

2.17 

1.046 

1.030 

0.3 

0.085 

2.74 

1.090 

1.054 

0.2 

0.144 

3.76 

1.185 

1.092 

0.1 

0.254 

6.18 

1.45 

1.166 

0.05 

0.351 

9.65 

I.89 

1.228 

0.025 

0.420 

11.42 

2.56 

1.270 

0 

0.510 

00 

00 

1.323 





MCA TM 1409 


3-5 


2. THE KHHMATIC CONDITIONS 


2.0 Statement of the Problem and Formulas 

2.0.0 General forms - apparent circular frequency 0 and phase 
velocity v at the front of the shock wave . - Let us consider in one of 
the regions Eq or E-j_ an acoustic field incident on S which we shall 

represent by the velocity potential 

<&(r,y,z,t) (2.0.1) 

with 4 and x having the subscripts (0 or l) corresponding to the 
region concerned. 

This field produces at the front S ' a phenomenon which with respect 
to the axes oXyz fixed to S will be governed by the form 

='®(at,y,z,t) (2.0.2) 

where a has the same subscript as • $. 

As a result of the fact that our boundary conditions are in point 
form, there must appear at S for every field associated ■with, the 
incident field by reflection or refraction a corresponding phenomenon 
whose spatial distribution and development with time are governed by the 
form t|t. This condition of kinematic order is independent of the dynamic 
conditions which derive from a particular choice of the equations con- 
necting the values of 5, s, u, 8, taken from both sides of S. 

Since we limi t ourselves to the case where the incident wave is 
plane, progressive, and uniform, we have 


$ = & “ ( xcoscp+ys ±m.cp ) /<£] 


(2.0.3) 


with 9 being the angle of incidence. 

The function t|t which represents the phenomenon at S then assumes 
the form 


i)r = e 


ifl(t-y/w) 


(2.0.4) 



l6 


MCA TM 1409 


with 


— = 1 - A cos cp = — — + .... ( J- ., . + . A l l£ (2.0.5) 

® 1 + T d 


w = fl _ 1 - A cos cp _ 1 A + (l + A)t^ 
c = cd sin cp sin 9 St ~ 


( 2 . 0 . 6 ) 


where 


T 



c, A, cd, cp, r should have the subscript of the region which contains 
the wave. 

2.0.1 Condition necessary for a plane, progressive, and uniform 
acoustic wave to be an incident wave . - Let us consider a mobile parti- 
cle M carried along by the waves and advancing on a sound ray. If M 
approaches the S plane with increasing t, the wave is incident. If 
it recedes from it, the wave is reflected or refracted. 

With respect to the axes fixed to the fluid in the region considered, 
the velocity at which M is displaced toward increasing x is c cos cp. 

In the region under subpressure Eq, the distance of M from S 
is, except for a constant, (cq cos cpQ - aQ)t. It decreases therefore 
constantly regardless of what cpQ may be since we have &q > Cq. Con- 
sequently, any plane, progressive, and uniform wave present in Eq is 
an incident wave. 

In the region under superpressure E-j_, the distance from M to S 
is, except for a constant, (a-j_ - c^_ cos cp-^t. 

Since a^ is less than c-]_> two cases are possible 


(a) 


a l " 


cp x < 0 where |cp^| < arc cos A-]_ (2.0. 7) 


C-j_ cos 


or 


rj< \J(1 - A 1 )/(l + A X ) 
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The wave is incident. 


(b) a^_ - C]_ cos cpj_ > 0 where arc cos A^ < jcpj < rt (2.0.8) 

or | T l| > \f0- “ A l)/( 1 + A l) 

The wave is produced by reflection or refraction. 

2.0.2 General relations .- The plane waves, whether acoustic or 
thermal, are governed by the form 

e ht-m(ox+Py) (2.0.9) 

with 


and 


m = ±h/c acoustic waves 

( 2 . 0 . 10 ) 

m = ± \j yh / k thermal waves 

On the S -plane which is displaced at the velocity a in the region 
considered, the exponent of (2.0.9) must assume the form i0( t - y/w); 
hence the relationships. 

i'ft = h - ama ifl/w = pm (2.0.11) 

that is 

a = (h - ifl)/am p = ift/wm (2.0.12) 

a 2 + p 2 = 1 gives 

(am/o) 2 + (a/w) 2 - (h/o - i) 2 = 0 (2.0.13) 

Replacing, in this equation, m by its value taken from one of the 
relationships (2.0.10), we obtain an equation which will determine h as 
a function of 0 and w. We then take m, a, and p from (2.0.10) and 
(2.0.12). 
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2.1 The Reflected or Refracted Thermal Wave 
Inserting into (2.0.13) the value m^ = 7 h/k gives 


(h/ft ) 2 - 2 (i + 7 a 2 / 2 kft)h/ft - (l - a^w 2 ) = 0 ( 2 . 1 . 1 ) 


whence 


^ = i + zaf 

ft 2 kft 


( 2 . 1 . 2 ) 


In order to he suitable for our purpose, the root - must have a 
negative real part: a wave which vanishes when t increases indefinitely. 

As we have seen, the dimensionless number kft/a 2 is very small. On 
the other hand, ( 2 . 0 . 6 ) gives ft/w = (cu/c) sin cp (where ca, c, and cp 
have the subscripts of the region which contains the incident wave). 

Consequently, kft/aw = (kw/ac) sin cp is the same order of magnitude 
as kft/a^, that is, very small. 

If we then take the approximate value of the radical of the formula 
( 2 . 1 . 2 ), we see that the suitable root is the one where this radical has 
the negative sign. 

In order to obtain an approximate value of h /ft, we must expand the 
radical up to the terms of the second order; the terms of the first 
order disappear and we obtain 


kft 2 a 2 + w 2 
7 a 2 w 2 


(2.1.3) 


One can see easily that 

kh _ 1 k% 2 (a 2 + w 2 ) 

c 2 7 c 2 a 2 w 2 

remains always very small. The approximation made in order to separate 
the thermal wave from the acoustic wave is therefore valid. 

On the other hand, (2.1.3) shows that h is real and very small. 
Consequently, the thermal wave is not perceptibly oscillating and its 
damping is slow. 
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From the value of h we derive 


m 


in \ja 2 + w 2 / 


aw 


( 2 . 1 . 1 )-) 


a , = cos cp + 

\] a 2 + w^ 


( 2 . 1 . 5 ) 


P = - f= ■= = sin cp t 

V a 2 + w 2 

The spatial distribution of the temperature excesses and of the 
condensations is therefore sinusoidal. The wave is oblique. It is not 
progressive and dies out on the spot. 

As a result, this wave can originate only behind the S-plane (the 
movement of which it cannot precede). Thus it can be present only in 
the region under superpressure E]_. 

The absence of any thermal wave in the Eq region justifies there- 
fore the method by which we have established the boundary condition 
( 1 . 1 . 20 ). 

The potential Q = e -m ( ax+ Py ) is then written 


Q = e 12 ( x A-y/v) (2.1.6) 

The quantities fixed to the wave then sire (see 1-2.13) 

H - " £ Q eht e t = -8 t ®t - 0 
and the velocity components 

"t - - 1 s t v t = 1 s t (2.1.7) 

These values will be useful for simplification of the dynamic 
conditions . 
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2.2. The Reflected Acoustical Wave 

In this case, the incident wave is, like the reflected wave, situ- 
ated in the region under superpressure E^. The relationships (2.0.5) 
and (2.0.6), connecting the kinematic characteristics cd and qp (or t) 
of the acoustic wave with the circular frequency 0 and the phase 
velocity w which can be observed at the front of the shock wave, are 
simultaneously valid for the incident wave (a>^> cp^_, tj\ and the 

reflected wave cp-j^, t-[) . 

Eliminating fi and w, we obtain 

T l T i - (1 " A x )/(l + A ± ) (2.2.1) 


sin cp 1 (l - AjJ^ + (l + Aj~)t^ 

“i sin (i - a x ) (i + a 1 )(i + tJ) 


(2.2.1) shows that to the incident wave 
there corresponds indeed a reflected wave Pr^ 



(2.2.2) gives the ratio of the circular frequencies (perceived by 
an observer carried along by the fluid) of the incident wave and of the 
reflected wave. 


When qp-j_ increases from zero (normal incidence) to arc cos A^_ 


(limiting angle of incidence), t-j_ increases from 
cd i/o>L increases from ^1 - A ^ j ^1 + A^ to 1. 


Thus we have 



and 


sin cpj^ 1 + A^ 

sin q>]_ < 1 - Aj_ 


Due to the fact that any plane, progressive, and uniform incident 
wave produces a reflected wave of the same type, Huyghens ’ construction 
is valid. 
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2.3. The Refracted Acoustic Wave 

2.3*0 Generalities - separation of the two types of refracted waves . - 
In order for refraction to occur, the incident wave must he s ituated in 
the region under subpressure Eq. On the other hand, any plane, pro- 
gressive, and uniform wave situated in Eq is an incident wave and will 
give rise to the production of refracted waves. 

The phase velocity w of the phenomenon produced on S by the 
incident wave is given by (2.0.6) 


= (Ap + - (Aq - X ) c 

2t 0 0 


(2.3*1) 


where 


Aq> 1 

When the angle of incidence cpQ varies from nt to +jc, Tq varies 

from -00 to +». w increases constantly and its variation is given by 
the following table: 



The same phase velocity w can therefore be produced by two waves; 
the incidence of the first lies between -:t and 0, that of the second 
between 0 and +n, with the first being larger in absolute value than 
arc cos 1/Aq, the second smaller than that value. The corresponding 

values of t, Tq, and r^, relative to these two waves, are connected 
by the relationship 
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T 0 T 0 = - (A 0 - 1 V(A 0 + 1) (2.3.2) 

Finally* (2.0.6) shows us that the circular frequencies of the two 
incident waves which produce on S the same phenomenon (equality of w 
and of fl) are in the ratio 

<* 6 /% = l sin V 8in *61 

Thus* the acoustic phenomenon observed on S* and all the more in 
the interior of the region* may be produced indifferently by one 

or the other* or any linear combination, of the two incident waves 
described above. Not a single observation made in E-, can serve to 
differentiate . 

The character of the refracted wave depends on the absolute value 
of the phase velocity w. If | w | is greater than ^jc^ - a^* the 

refracted wave is plane* progressive, and uniform. If |w| is less 
than one cannot make any plane* progressive* and uniform 

wave in E-j_ correspond to it. The refracted wave then assumes an 
exponential structure and we shall call it accompanying wave. 

The angles of incidence cp^ corresponding to the value of w* at 
which one passes from one type of refracted wave to the other* are given 

"by 




2 , 


®0 + c 0 



(2-3.3) 


When |tq| is outside the interval defined by the above boundaries* 

the refracted wave is plane* progressive* and uniform. In the opposite 
case we obtain the accompanying wave. 

Expressed as functions of the theoretical relationships which 
govern the shock wave* the two limiting values of r Q are written 



ti(l - £) ±\jl + HI t \f(n + 1)(1 +~T) 

n + i \JiTTT +-\](V + i)i 


(2.3.4) 
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The values of the corresponding limiting angles are given in 
figure 1. 

2.5.1 Plane, progressive, and uniform refracted wave; w^ > c|[ - a^ . 

The equation (2.0.6) applied to the two regions Eq and E]_ which con- 
tain the incident wave end the refracted wave, respectively, gives the 
relationship which connects and 


2 2 
C 1 ” a l + ( c i + a l) T l c 0 ” ®o + ( c 0 + a o) T 0 


(2.5.5) 


Since Tq is governed "by the incidence 9 q, is determined hy 

an equation of the second degree of which only the root, which is larger 
than \J( c j_ ~ a l)/(c]_ + a l), corresponds to the refracted wave. The 
other root gives in the incident wave which could, hy reflection, 

produce the refracted wave. 

Thus we have 


T n = 


a 0 + ( c 0 + a ' 


o) T 0 + \l[ c 0 - *0 + ( c 0 + a o) T 0_ 2 " K C 1 " a l) T 0 


2( c i + a i) T 0 


(2.5.6) 

which we may write , using the dimensionless functions Aq, A-^, r. 


1 - 


T_ - 


Aq + (! + Aq) Tq +^[( 1 - Aq) + (1 + Aq)t§] 2 - ^ (l - a|)t§ 


1 + At 
2 1 T- 

r 0 


( 2 . 5 . 6 *) 

The value of the circular frequency of the refracted wave, 

perceived hy an observer who is carried along with the fluid in the 
region E^, is then given hy (2.0.6) whence we derive 
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01^ sin <Pq 
^ “ r sin cp i 


(2.5*7) 


For an observer fixed to the axes, carried along with the fluid 
situated in Eq, there occurs a shift in frequency and the circular 

frequency perceived will be cjd£. 




CD. 


U 


1 = 0i l[ 1 + ^ COS 'Plj 

sin Oi, 


( 2 . 5 . 8 ) 


= 1 + — COS 9 - 


. , Q. J2 

1/ sin 9, T 


2.5.2 The accompanying acoustic wave w^ < c 2 - ai^.- We shall show 

that the refracted wave may then be represented by a velocity potential 
of the form 


# = e (lxo 1 -8 1 )t+[(iTi+C)x+(iTi , +^')^c 1 


must satisfy the general equation of sound 
which involves the relationships 



C 2 - n 2 + £’ 2 - 11 « 2 = s| - cof ^ , -Bf^ (2.3.10) 

On the other hand, on the front S, that is, for x = a-jt, the 
exponent of J>]_ must take the form ift( t - y/w) whence 

a»l - 81 + ^ (iT) + 5) - ift ^ - 0^ - - -1 § (2.5.11) 

From (2.5-10) and (2.5. ll) we then take 


v° - VC 1 - A !) 


(2.5.12) 


8 

a 




(2.5.15) 
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Ti/n = -a-l^i - a|) = -a-jc 


(2.5.14) 


i = ii 

£2 £1 


(2.5.15) 


t[ y& = -c^w 


(2.5.16) 


S’ = 0 


(2.5.17) 


The equation (2.5-15) shows that this type of wave exists only if 

the phase velocity w is smaller than ■= 1 \4 - 4 - nFT 7 a^j conse- 
quently, this wave appears only in the case studied in the present 
section. On the other hand, this condition of existence implies that 
we have > a-j_. No wave of this type can therefore appear in the Eq 

region where Cq is less than aQ. 

Before describing the accompanying wave, we must verify that the 
condition * a « 1 is satisfied, that is, that we have always ks/ cfj- « 1. 
Replacing fi in (2.5-15) "by its value taken from (2.0.6), we obtain 


kS , r x ) 1 - *1 - (* 7 °!^ 

-P 1 - .P 


sin cp Q M 
c 0 


(2.5.18) 


ktn/cQ is very small) A^_ and r lie between 0 and 1) thus we 
have a very small ks/c^. 

In order to make the characteristics of the wave evident, we shall 
take the axes oXY whose origin o is carried along with the S -plane 
and is displaced toward increasing y at the velocity w. We therefore 
set 


X = x^ - a^t 


Y = y - wt 


(2.5.19) 
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The potential then takes the form 

- £ (l ^W3/=l (2.5.2C 

which Is independent of t. 

The wave is therefore motionless with respect to the axes oXY. 

(;x/ c i 

Its amplitude decreases as e if one moves away from the shock 

wave (toward the negative x). This wave constitutes therefore a 
phenomenon which runs behind the shock wave while sliding at the same 
time transversely at the velocity w. These characteristics justify 
the name we gave it: accompanying wave. 

Let us define its structure in detail. Toward increasing Y the 
phenomenon is spatially periodic and its wave length is 


Ay = - 2jtw/ft 


(2.3.21) 


which was, besides, evident in view of its origin. Toward decreasing’ 
X the phenomenon is spatially of a damped pseudoperiodic nature. Its 
wave length is 


Ay = 2rtCy/T] = 2rt (c^ — a lj/ a l^ (2.3*22) 

The planes of equal phase, or pseudo-wave planes are given by: 
rjX + q'Y = const. They form, therefore, with the S-plane of the shock 
wave the angle X given by 

tan X = = “( c i “ a ?)/aiv (2.3.23) 

This formula permits a very simple geometric construction Of that angle. 
The rate of damping of the wave toward, the negative X is 


x 

C 1 



(2.3.24) 
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It is zero for the limiting values w = ±u c^ - a^. It increases when 
w decreases. v 

For w = 0, we have on the one hand £1 = 0, on the other 
fl/w = 0 )q sin cPq/Cq and cos cpQ = Cq^Sq 
^ l/ C l ^ en assumes the value 


= c i J4 - H 

W/v -0 ^ . a 2 =0 =0 


which may be expressed as a function of the variable | which defines 
the amplitude of the shock wave* and of the wave length 7 \q = 2jtc^Dq 
of the incident wave in the region Eq 



2tc f ii~+ 1 
+ 6 


(2.5.25) 


In order to calculate the characteristics of the accompanying wave 
from those of the incident wave, it suffices to replace ft and w in 
the relations (2.3.12) to ( 2 . 3 • l6 ) hy their values taken from (2.0.5) 
and (2.0.6), namely 


n - ^(co - ®o COB < Po)/ c o 
w ■ ( c 0 ' ®o cos 'Po) / sin ^0 

J 

The formulas obtained are complicated. It is unnecessary to write them 
here explicitly. 


2.4. Geometric Constructions 

2.4.0 Huyghens 1 construction (fig. 2) .- With the plane of the 
figure being motionless with respect to the region under subpressure 
Eq, S and S' are the positions of . the shock-wave front at the times 

t » 0 and t = 1. The distance SS* is slq. . . 
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Let us consider in the region Eq, to the right of S, a plane, 
progressive, and uniform wave which intersects S at the point 

at the time t = 0. At the time t = 1, this wave will intersect S' 
at a point D and its potential prolongation will be tangent to a 
circle Cq, with the radius cq and the center Cfo, entirely situated 

to the left of S’ since c q < a Q* 

Assume B to be the foot - of a line through 0 q perpendicular to 

S'. BD will be equal to the phase velocity w observed on the shock- 
wave front . 

We can always make two tangents to the circle Cq pass through D. 
The figure therefore confirms that in the region Eq two waves Iq, 1^ 

exist to which the same phase velocity w corresponds on S. These two 
waves are incident because the point of contact of the tangents drawn 
from D to Cq is to the left of S: that is the region where Eq is 
virtual. 

Let us now consider the region E^. The fluid flows there toward 

the right, with the velocity U. A- plane wave passing at the time 
t = 0 through the point 0 0 , will, at the time t = 1, be tangent to. 

a circle with the radius c^, whose center 0^ is at the distance 

OqO-^ = U from Oq. 

O-jjB is therefore the velocity a^ of the shock-wave front with 
respect to the fluid filling up the region Ej_* Since a-j_ is lesB than 
c-p the circle intersects the straight line S' at two points, F 

and F 1 . 


The incident waves are those whose point of contact with C-j_ 

is to the right of S' (there where E-j^ is virtual); the reflected or 
refracted waves touch to the left of S'. 

Let us now assume D to be a point on S' situated outside the 
segment FF'. From this point we may draw two tangents to C^, one of 

which corresponds to an incident wave while the other corresponds 

to a reflected or refracted wave R^. 

Thus, if we assume a wave R^ in E-^, this wave determines the 
point D, and we see that it can have been produced either by reflection 
of the wave 1^ (situated in E^) or by refraction of one or the other 
of the waves Iq and 1^ situated in Eq. 
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Finally, it appears that when the point D lies "between F and 
F’ , there correspond to it two more possible incident waves, Iq and 

Iq in Eq, but none in since D lies inside the circle C-j_. 

The cross-hatched sectors of the figure show the angles of incidence 
cpQ for which this condition is realized. Huyghens' construction thus 

does not permit us to know what is happening in the region E]_. 

2.4.1 Construction of the accompanying wave (fig. 2) .- Formula (2.3.23) 
gives the angle X which is formed by the pseudo-wave plane (plane of 
equal phase) of the accompanying wave and the S-plane. It lends Itself 
to a simple construction. 

Let us draw the tangent from the point 

It intersects the straight line O^B at G. 

Since O^B is equal to a-j_, we shall have 

■ 4 

20 “ ■§£ " a l 

Since, on the other hand, D^, situated between F and F* , is the 

point where the incident wave intersects S’ at the time t = 1, we 
have BD-j_ = w. Consequently 


F’ to the circle Cq_* 

0-lG = O^/o-lB = c^/a-L- 


BG 

BD X 


c l/ a l " a l 
w 


tan X 


Thus, the angle BD^G is the desired angle X. The prolongation 

of the straight line GD-j_, to the left of S’, gives the pseudo-wave 
plane of the accompanying wave. 

The construction thus shows very clearly the variation in direction 
of this pseudo-wave plane as a function of w, that is, of the position 
of the point D-^. 

Figure 2 shows these different constructions for a whole series of 
incidences. The arrows placed on the traces of the various waves indi- 
cate the direction of their propagation. The sketch corresponds to the 
case where the shock wave involves an abrupt doubling of the pressure, 
that is, where one has 


Po/Pi = 6 = 0.5 



50 


NACA TM 1409 


2.5 Use of the Acoustic Phenomena for the Study 
of the Shock Waves 


The kinematic conditions are expressed by formulas into which 
enter only velocity ratios: sonic velocity, velocity of the wave front, 

flow velocity. They do not stipulate the validity of the theoretical 
relationships which express these velocities as functions of the shock 
amplitude . 

Moreover, in the case where the shock-wave front is preceded or 
followed by compression or expansion phenomena which entail variations 
in temperature and in the flow velocity of the fluid, the acoustic waves 
undergo continuous refractions. These phenomena, incidentally, are well 
known. 

Consequently, the exper im ental study of the kinematics of the 
acoustic waves reflected or refracted by a shock-wave front furnishes 
a research method which is probably not without interest for aerody- 
namicists . For undertaking it, we may consider photographic, even 
motion-picture methods which show the form and the deformations of wave 
fronts, or microphonic methods which detect the variations of frequency. 

One will note that a large amount of information can be collected 
solely by the study of the reflection of waves, that is, by means used 
exclusively in the region under superpressure E^. « 


3. DYNAMIC CONDITIONS 
3-0. Elimination of Thermal Waves 


Let us start from the relationships established in section 1.1.2, 
namely 


01 - 9 0 = 2| £( 5 ^ - m 0 ) (1.1.17) 

8 1 - s o " £ |- (“1 - E o) 


(1.1.18) 
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c l u l 




( 1 . 1 . 21 ) 


We have seen that in our problem the region under subpre&sure E q 
can contain only incident waves, representable by the velocity potential 


1^0 t-(xocoscp 0 +y 0 sin<p 0 ) /c Q 


( 5 . 0 . 1 ) 


Thus we have 


s o - - 1 ( a ^/ c o 2 )*o 

“0 = 7*0 ~ -ir(a^/c§)« 0 . 

e 0 “'(y - l) s 0 * - ±( <7 - l)(^o/ c o)®0 

Uq = -i(^o/ c o)cos cp 0 $ Q = c Q cos cp 0 s 


( 5 . 0 . 2 ) 


In the region E^_ there is a superposition of an acoustic phe- 
nomenon (letters with the subscript la) and a thermal phenomenon 
(subscript It). 


The formulas (2.1.7) give us for the thermal phenomenon 


“it " 0 e it = " s lt u lt = ~ i ( kn /y a l) s lt (5-0.3) 

In the acoustic phenomenon we have 


“la ■ 7 s la (5-0. ^) 

and in all 

“1 " “It + “la " “la 9 1 = e lt + 9 la 


S 1 " s lt + 8 la U 1 ■ u lt + u la 
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The equation ( 1.1,18) thus yields 


s lt = "7 J (“i - & 0 ) (5-0.5) 

where J is the function defined "by the formula (l.l.ll). 

This relationship will give us the thermal wave as soon as we have 
calculated the acoustic wave. 

Introducing this value s 1+ into the third relation (5-0.3 )t we 
have 1 


u lt = - ® 0 ) 


( 5 - 0 . 6 ) 


This we insert into (1.1.21) which becomes 


c l u la 


- 0 - i 


kfl 

a l c l 


7J, 


2 - 

°l“l _ 




- H - i 


kfl 
an Ci 


I'-O 


7Ji 


2 - 

c cPo 


(5-0.7) 


We know that the numbers kfl/a^c^ and kH / a i c o are very small. 

On the other hand, table II in section 1 shows us that the function J 
is always smaller than the functions G and H. Consequently, we can 
neglect in equation (3-0.7) the Imaginary terms, that is, the terms 
stemming from the thermal wave. 


Thus there remains a condition which involves only the acoustic 
waves present on both sides of the 'shock-wave front S 


c l u la 


G 


c l cu l = 1 

7 r 




( 5 - 0 . 8 ) 


3-1. The Reflected Wave 

We have to consider in the region two waves, one incident, 

the other reflected, represented by their velocity potentials 
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4l - 

*. , Bie 1 ®Lf-K“ 0 'l , i + yi Bin -’’i)/ c ll (3.1.2) 

The Eq region, on the other hand., does not contain any wave. The 
second term of the condition (3*0.8) is therefore zero. 

Assuming the kinematic conditions satisfied (3-0.8), we then have 


Bjpi (cos cp]_ - g) + Bjcn^cos q>j_ - G) = 0 (3-1-3) 


The amplitudes of the waves are proportional to the factors Bto. 
The ratio of the amplitude of the reflected wave and of the amplitude 
of the incident wave, that is, the amplitude-reflection coefficient, R, 
then is 


R = - (cos cp 1 - G) /(cos cp£ - g) 


( 3 - 1 -^) 


Let us refer to the kinematic relationship (2.2.1) which connects 
the tangents t-j_, of the half angles cp^, i we can eliminate 

cp-J_, and we obtain 


R = 




+ tS 


1 - 0 _ J2. 
.1 + G 


1 + 


'1 ~ AjY _ 1 - G t 2 
,1 + An / 1 + 0 1 


( 3 - 1 - 5 ) 


with < (l - A^)/(l + A^) (incident wave). 


For cp]_ = 0, t-j_ is equal to 0 (normal incidence), Rq is equal 

to (l - G)/( 1 + G), which is negative because G is greater than 1. 

(See table II, section 1.) 

For <p^ = arc cos A^, t-j_ = ± \j (l - A^)/(l + A-^ (limiting incidence), 
R lim is e< l ua l to "1- 
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Since we have in this case = <f>n. 


= 6D i-’ value signifies 


that in the region where the two waves are superimposed on one another 
(and which, incidentally, is vanishing) their sum is zero. In other 
words, the phenomenon involves only an incident acoustic wave which 
stops at S, which is actually the case. 


Thus, the coefficient of amplitude reflection is constantly nega- 
tive and, in absolute value, increases from (G-l)/(G+l)to 1 when 

q> 2 . increases from 0 to arc cos A^. 


Thus it will he possible, in practice, to observe the reflected 
wave at incidences close to the normal only if the shock wave is intense 
(G clearly larger than l). 


3.2. The Progressive Refracted Wave 

Let us represent the incident wave, situated in Eq, by its velocity 
potential 


= B 0 e 


^Op"( X 0 COSCp 0 +y 0 8lnC Po)/ c oH 


and the refracted wave, situated in Ep by 

ico.jjt- (x 1 coscp 1 +y 1 sin«P 1 ) j c^j 


= B-ie 


Starting from the instant where the kinematic conditions are 
satisfied, the exponents of the exponential terms are the same on S. 
The dynamic condition (3.0.8) then gives 


IB-[C0i(c o s q>! - G) = B o“o( cos cp Q - h) 


( 3 - 2 . 1 ) 


We shall compare the amplitudes of the acoustic pressures in the 
two -waves, that is, we shall form the ratio 


R = PA/Po^ = PpB^/PqBcPo 
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the value of which is given "by (3*2.1) : 



H - cos cp Q 
G - cos <P]_ 


where R is equal to Pq/P-j_ and 


( 3 * 2 . 2 ) 


jl = IIjl±J = r 

Rr Mid + hi) i 

(See table II, section 1.) 

In order to calculate R, one must refer to the kinematic conditions 
which give <p^ as a function of 9 q* 

We shall give below only the values R as functions of | for 
<Pq = 0 (incident wave parallel to the shock wave and receding before 

the latter) and <Pq = it (incident wave parallel to the shock wave and 

coming to meet it). 

We have in the two cases cp^ = jt whence 


cp = o R = — 5 -Z-J: 

^0 K 0 Rr G + x 


p = 1_ H + 1 

^ Rr g + i 


cp 0 = rt 
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i 

R 

for 

O 

II 

s? 

q>0 - it 

1 

0 

1 

*9 

.0002 

1.094 

.8 

.0010 

1.22. 

*7 

*0035 

I.36 

.6 

.0092 

1*55 

*3 

.022 

1.82 

.4 

.049 

2.19 

*3 

.120 

2.79 

.2 

*332 

3*92 

.1 

1.29 

7.00 

*05 

3*85 

12.5 

.025 

7.85 

17*9 

0 

00 

00 ■ 


3*3* The Accompanying Wave 

With the incident wave being represented in the foregoing manner, 
the accompanying wave is derived from, the potential 

*1 “ B 1 i i£Dl “ 8 l) t+ P ±n+? )Xl+±11 /°1 

When the kinematic conditions are satisfied, the boundary condition 
(3*0.8) gives 


H - cos <p. / . u \ 

ifflQBo 53— 2 - (*“l - 8 l)Bl (0 -' jg! ij (3.3.1) 

Ihe ratio R which interests us most is the ratio of the acoustic 
pressure immediately behind the shock wave p^ = - &j)® 1 an< ^ 

of the acoustic pressure in the incident wave Pq = -Ip^qBq. 

Taking the relations (2.3*12 to 2.3*17) into account, we obtain 

n _ cos 9 q - H i - 

- m 1(*1 - G) - (^^(l/Ai - G) 


(3.3.2) 
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and from the fact that 
obtain 


5i/o i is equal to A]_ ^ ^1 - A^j ( c i / w ) 2 " 1* 


we 


R = 


H - cos cp c 
IR 


I 


Af - W<i 


(l - A x 0 ) 2 + (G 2 - l)wf 


where 



r 1 - A o °°° ’Pq 

sin cPq 


For illustrative purposes we have performed the calculation of | R| 
in t he extreme case where w is equal to 0, that is, cos <Pq is equal 
to 1 /Aq. We have then 


I R iv=0 


i » - y*o 

TR ljA 1 - G 


For cpQ s arc cos 1 /Aq we obtain 

5 = 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.05 0.025 0 
I r | = l 1.05 1.15 1.31 1.46 1.71 2.0 6 2.66 3.86 7.45 l4.6 22.8 eo 

We see that, as in the ‘previous case, the acoustic pressures may 
be increased by a significant factor. 


4. SUMMARY OF RESUME 


The presence of acoustic waves in one or the other of the regions 
which separate the shock-wave front is manifested in the following 
phenomena: 

(a) In the region under subpressure Eq in which the shock wave 
progresses at a supersonic velocity, no type of reflected or refracted 
wave is produced. 
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(b) In the region under superpressure E-j_ the encounter of an 
incident acoustic wave (coming from the Eq region or the E^ region) 

with a shock wave causes the formation of two waves of different type. 

One is an acoustic wave (refracted or reflected depending on the origin 
of the incident wave); the other is a thermal wave. 

(c) The thermal wave is not progressive in the first approximation 
(that is, when the incident acoustic waves are of an ordinary amplitude). 

It is therefore motionless with respect to the gas contained in the 
region under superpressure and carried along by this gas in its motion 
of overall flow. It then moves away from the shock-wave front. This 
wave involves spatial periodic variations in temperature and in density 
but no variation in pressure. It is damped slowly and aperiodically. 

(d) When the incident acoustic wave is situated in the region under 
superpressure, in that same region a reflected acoustic wave, of the 
same type as the incident wave, is produced; that is, in the case 
studied here, a plane, progressive, and uniform wave. 

Huygens’ construction permits obtaining its kinematic characteristics: 
direction of propagation, and wave length. 

The amplitude of the pressure variations in the reflected wave is 
always smaller than the amplitude of these variations in the incident 
wave. 


(e) When the incident acoustic wave is s ituated in the region under 
subpressure, a refracted acoustic wave is produced in the region under 
superpressure; the characteristics of the latter wave depend on the phase 
velocity w of the phenomenon produced on the shock front by the inci- 
dent wave. The main characteristics of the refracted wave are the 
following: 

(a) In the region under subpressure Eq there exist always 

two incidences of plane waves which produce on the shock front the 
same phase velocity w and, consequently, in the region under 
superpressure the same refracted wave (except for the ampli- 

tude ratio). 

The observation of the refracted wave in the region under 
superpressure E]_ therefore does not permit distinguishing 

whether this wave has been, produced by one or by the other of 
the two possible incident waves, or by a combination of both. 

( p) If the phase velocity w of the disturbance produced on 
the shock wave by the incident acoustic wave is larger than a 
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certain limiting value Wq (which depends on the intensity of the 
shock), the refracted wave is of the same type aB the incident 
wave: plane, progressive, and uniform. Huyghens 1 construction 

permits obtaining its kinematic characteristics. 

The amplitude of the pressure variations in the refracted 
wave may he, depending on the case (angle of incidence, shock 
amplitude) smaller than, equal to, or larger than the amplitude of 
these variations in the incident wave. All things being equal, 
this amplitude increases with the shock amplitude. 

The increase in amplitude may become very important in the 
case of intense shocks. 

(y) If the phase velocity is smaller than the limiting velocity 
Wq, which corresponds to incidences situated in a sector the 

boundaries of which are functions of the shock amplitude, the 
refracted wave assumes the form of an "accompanying wave." This 
wave remains attached to the shock wave which it accompanies while 
shifting parallel to that wave at the velocity w. 

The wave planes (planes of equal phase but not-equal intensity) 
are oblique with respect to the shock. Their direction may be 
obtained by a simple geometric construction. 

The wave amplitude decreases exponentially as one moves away 
from the shock-wave front. Immediately behind the shock wave, the 
amplitude of the pressure variations in the accompanying wave 
increases with increasing shock intensity. It may be larger than 
the amplitude of the pressure variations in the incident wave. 

(•8) In summary, the presence of acoustic waves in the region 
under subpressure, in which the shock wave progresses, may mani- 
fest itself behind the shock wave, on the one hand, by a field of 
progressive acoustic waves which expand gradually in the entire 
region under superpressurej on the other hand, by the appearance 
of accompanying waves, that is, of oscillating disturbances which 
the shock wave carries along behind it in its progress. 

Both may, in certain cases, involve oscillations of pressure 
very much larger than those that existed before the arrival of 
the shock wave. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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